Laser cooling of a trapped two-component Fermi gas 
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We study the collective Raman cooling of a trapped two-component Fermi gas using quantum 
master equation in the festina lente regime, where the heating due to photon reabsorption can be 
neglected. The Monte Carlo simulations show, that 3D temperatures of the order of 0.008 Tp can 
be achieved. We analyze the heating related to background losses, and show that our laser-cooling 
scheme can maintain the temperature of the gas without significant additional losses. 
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The realization of Bose-Einstein condensation (BEC) 
PJ and degenerated Fermi gases Q has outbursted the 
interest in the physics of ultracold bosonic and fermionic 
gases. A Fermi gas with attractive interactions under- 
goes for temperatures T below a critical one, T c , a transi- 
tion into the superfluid Bardeen-Cooper-Schrieffer (BCS) 
phase Q . The accomplishment of BCS and its possible 
detection have been considered in detail (see e.g. Refs. 
0)- Unfortunately, evaporative cooling currently em- 
ployed in experiments is based on collisional processes, 
and has not yet allowed to reach T < T c , since T c is 
much smaller than the Fermi temperature Tp, and for 
T < Tp the collisions are strongly suppressed due to 
Pauli blocking Recently, however, it has been pro- 
posed that T c may be significantly increased by employ- 
ing Feshbach resonances 0, or optical lattices Q- This 
opens the possibility to achieve BCS at a temperature 
where the collisions are still efficient enough. 

In this Rapid Communication we study the laser cool- 
ing of a two-component Fermi gas in the festina lente 
(FL) regime @, where the spontaneous emission rate 7 
is smaller than the trap frequency uj. In this regime the 
heating due to photon reabsorption is prevented, as re- 
cently observed for tightly bound atoms in optical lat- 
tices 0. The laser cooling in the FL regime has been 
already predicted to work for bosons and for polar- 
ized fermions ^]J. By using Monte Carlo simulations, 
we show here that laser cooling allows for bringing a two 
component Fermi system to T -C Tp. 

The laser cooling of fermions toward T <C Tp is obsta- 
cled by several problems. One of them is the inhibition 
of spontaneous emission, which results in the decrease 
of the cooling efficiency. In Ref. [llj we have predicted 
that the inhibition problem can be overcome by either 
dynamically adjusting the spontaneous emission rate in 
a Raman cooling process, or by employing specially de- 
signed anharmonic traps. In this paper we apply the for- 
mer solution. Another obstacle is related to the inelastic 
losses that create holes deeply in the Fermi sea, produc- 
ing a significant heating. We show that background col- 
lisions do not affect significantly the laser cooling. In 
fact, our cooling method can be employed to maintain a 
two-component Fermi gas at a fixed T, for a relatively 



long time in a trap. Finally, we derive an analytic for- 
mula that describes the heating of the trapped gas due 
to background collisions. The heating is shown to be 
smaller than in the homogeneous case. 

We consider fermionic atoms with an accessible elec- 
tronic three-level A scheme, containing states \g), \e) and 
\r). The ground state \g) is coupled via a Raman transi- 
tion to the metastable state |e). The latter state is also 
coupled by an optical transition to the upper state \r), 
from which atoms rapidly decay into \g). The adiabatic 
elimination of |r), leads to an effective two-level system, 
characterized by tunable parameters: Rabi frequency f2, 
and spontaneous emission rate 7. The value of 7 can be 
controlled by modifying the coupling from |e) to \r). 

The atoms are placed in a dipole trap characterized by 
a Lamb-Dicke parameter r\ — 2na/\, with a = y/K/2mw 
being the size of the ground state of the trap, and A 
the laser wavelength. A dipole trap was recently used 
for the all-optical production of a degenerate gas of two 
Li species and in current experiments in Mg 
We consider a spherically symmetric trap with incomen- 
surable frequencies oj 9 and uj e , for the ground and the 
excited state respectively. The latter assumption simpli- 
fies the dynamics of the spontaneous emission processes 
in the FL limit. The cooling process consists of sequences 
of Raman pulses of frequencies, adjusted in such a way 
that they induce the transition of atoms to the lower 
motional states of the trap. We assume that the power 
of cooling lasers is sufficiently weak, and therefore dur- 
ing each pulse no significant population in |e) is present. 
This allows to eliminate adiabatically the level |e), and 
to consider only the density matrix p(t) for the atoms in 
\g), and being diagonal in the Fock representation corre- 
sponding to the bare trap levels. 

We assume that the laser acts only on one component, 
whereas the other one is cooled sympathetically, which is 
sufficient to reach T <C Tp. Using the standard theory 
of quantum-stochastic processes [3 0] we derive the 
quantum master equation (ME) for the density matrix 
p(t) in a similar way as that for bosons [To| : 

p(t) = £ p + Cip + £ 2 p, (1) 

where C P = -iH cS p(t) + ip{t)Hl s + Jp(t), dp = 
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-i[Hl m , p(t)], C 2 p = -i[H con , p(t)]. The first term in 
the ME (JIJ describes the evolution of the atoms in the 
vacuum of the electromagnetic field. The effective non- 
Hermitian Hamiltonian H c s is defined as 

m 1 

+ ^2u b n bl 1 b n -i^J2£, lm elg m gl n e h (2) 

n l,m 

and the so-called jump super-operator J is given by 

J P(t) = 6m# m eip(i)e}# m . (3) 

l.m 

Here g m (ei) is the annihilation operator of atoms of 
the first component, in the ground (excited) state, and 
in the trap level m (1). The annihilation operator of 
atoms of the second component, and in the trap level 
n, is denoted by b n . These operators fulfill the stan- 
dard fermionic anticommutation relations: {g m ,gn} = 
{e m , e+} = {& m ,fyU = <5m,n- In Eq. © 7 denotes 
the single-atom effective spontaneous emission rate, wj^, 
w m , cj m are the energies of state m of the trap: \g), 
|e), \b) respectively, and 5 is the laser detuning from the 
atomic transition. The coefficients £ i m are defined as 
follows £i m = f^dcj) Jq d9 sin6> W{9, (/>)\r]i m (k)\ 2 , where 
W(9, is the fluorescence dipole pattern, and r)i m (k) = 
(e, l|e zfc ' r |<7, m) the Frank-Condon factors. 

The interaction of the laser light with the atoms is 
governed by the Hamiltonian i?i as : 

l.m 

where fl is the Rabi frequency and &l is the wavevector 
of the cooling laser. Binary collisions are described by 

#coll= Y Um.n.q.pgLgnb^bp, (5) 
m.n,q,p 

where due to Fermi statistics only collisions between 
different species are allowed. We neglect the col- 
lisions with the atoms in |e), since only a small 
fraction of atoms is excited in each pulse. The 
collisional amplitudes are defined as ?7 m ni q iP = 
(Awha sc /m) J R3 d 3 £(/> m (x)^ n (x)/3*(x)/3 p (x), where (j) n 
(Pn) denotes the wavefunction of the state n in the \g), 
trap, and a sc is the scattering length. 
We require the FL condition 7 < u> at least in the ini- 
tial phase of the cooling process, and also assume that 
f2 < 7 and fi 2 /7 <C u>, which allows for the adiabatic 
elimination of |e), performed by means of projection op- 
erator techniques [14| . For weak interactions between the 
two species, the dynamics due to the collisions and the 
dynamics due to the laser-cooling are independent |16| . 
The laser cooling may be described by the ME without 
the collisional part, whereas the collisional part takes the 



form of a quantum Boltzmann master equation (QBME) 
[13 . The probability of a laser-induced transition from 
state m to state n of the trap \g) is given by 

7 

x \p 7 2 6n|??lm(fcL)| 2 

^ 4[<5 - ( W f - uL)] 2 + l 2 R m i 2 ' 1 ' 

where R m \ = ^ n ^in(l — Nn + S n ,m) are the factors mod- 
ifying the spontaneous emission rate in a Fermi gas, and 
and 7V m denotes the number of atoms occupying the 
state m of the trap \g) and |6), respectively. In the regime 
of quantum degeneracy, coefficients i? m i vanish, inhibit- 
ing the spontaneous emission and forcing the atoms to 
remain excited for a long time. This prolongs the cool- 
ing process, and can result in excited-ground collisions, 
which lead to heating and losses. In addition, the adia- 
batic elimination used in the derivation of Eq. © ceases 
to be valid. The cooling efficiency is also decreased due 
to the fermionic inhibition factor (1— JV n ) in the numera- 
tor of © ■ The negative influence of the statistics can be 
overcome by dynamically increasing 7 during the Raman 
cooling in order to avoid the inhibition effects, but 
still remaining in the FL regime: 7-R m i < Still, a small 
fraction of the atoms will remain in the excited state af- 
ter the cooling pulses, and has to be removed from the 
trap in order to avoid inelastic collisions. The latter aim 
can be achieved by optically pumping the excited atoms 
to a third non-trapped level. 

The probability of a collision between two fermions of 
different species, from the states n and p to the states m 
and q, respectively, is given by 

X \U m ,n,q,p\ 2 ^E n +E p ,E m +E^- (7) 

The fermionic inhibition factors (1 — N^) also slow the 
collisional processes. 

In our simulations we assumed the atoms as con- 
fined in an isotropic, 3D harmonic trap with frequency 
u> = oj 9 = u> e = oj b . Due to the limitations when sim- 
ulating relatively large systems (N ~ 10 4 ), we employ 
ergodic approximation, i.e. we assume that the popula- 
tions of the states with the same energy are equal. This 
approximation relies on the fact that thermalization in- 
side the same energy shell is much faster than between 
different energy shells. In addition we assume that the 
collisional processes is much faster than the laser cool- 
ing. The detailed calculation of the transition probabil- 
ities within ergodic approximation is presented in Ref. 
|18|. Due to numerical limitations we assumed rj = 2, 
which could be e.g. the case of potassium atoms in a 
dipole trap with lo = 2n x 2.4 kHz (employed in Li ex- 
periments [l2|l. and a laser wavelength A ~ 720 nm. 
We consider a scattering length for the interactions be- 
tween the two species a sc = 157ao, where ao is the Bohr 
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FIG. 1: Time-dependence of the temperatures of the two 
components (darker and lighter curves) during the cooling 
by controlling the effective spontaneous emission rate. Inset: 
time dependence of the number of atoms. The losses are due 
to background collisions and the removal of long-living excited 
atoms. For details, see text. 
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FIG. 2: Evolution of a degenerate two-component Fermi gas 
after the cooling process, in the presence of background colli- 
sions with a rate 7b g = 350 Hz. The temperatures of the two 
components are plotted (darker and lighter curves). 



radius. This value corresponds to the interactions be- 
tween]^ 1 = 9/2, m F = 9/2) and \F = 9/2, m F = 7/2) of 
40 K |2l|. The trap is assumed to have the same depth 
for both species, and contains 81 energy levels (91881 
states). We assume both components to have equal ini- 
tial number of atoms, N — 10660, corresponding to a 
Fermi energy Ep — i%hio. For this relatively large num- 
ber of atoms, the dynamics generated by the collisional 
part of the ME equation leads to an equilibrium distri- 
bution, which agrees very well with the one calculated 
from the grand canonical ensemble 01- Hence we start 
the simulations from a thermal distribution. 

In a Fermi gas the main loss sources are provided by 
background collisions, resulting from non ideal vacuum 
conditions, and photoassociation. The photoassociation 
losses (when the laser is tuned between molecular reso- 
nances) are typically of the order of 10~ 14 cm 3 /s for laser 
intensities of lmW/cm 2 [2(j. In our case the applied 
laser intensities are typically 1000 times smaller and we 
estimate that for N = 10660 atoms, the atomic density 
is smaller than 3.5 x 10 14 cm -3 . Therefore the photoas- 
sociation losses can be safely neglected. On the contrary, 
the background collision apart from decreasing the num- 
ber of atoms, generate holes deep within the Fermi sea. 



Those holes, after subsequent thermalization, may lead 
to a significant heating [lj| • We have assumed that back- 
ground losses depopulate each state of the trap with the 
rate 7b g , which is independent of the energy of the state 

= -7b g ^, (8) 

where j — 1,2 enumerates the components. In the sim- 
ulations we assume 7b g = 1/350 Hz 0. We include also 
in the simulations the losses produced by the removal of 
long- living excited atoms at the end of each cooling pulse. 

Fig. n shows the evolution of the temperature of 
the laser-cooled two-component Fermi gas. Initially 
we consider Tq — Tp. The cooling process was di- 
vided into four stages, each one consisting of a se- 
quence of two Raman pulses. The employed pulses 
are characterized by the following parameters: detuning 
S/u> = {(-11, -12), (-16, -17), (-19, -20),(-25, -26)} 
respectively, Rabi frequency = {(0.113,0.113), 

(0.008,0.012), (0.0025, 0.004),(0.003, 0.01)} respectively, 
and length At/iv' 1 = {(250,250), (2000,2000), 
(4000, 4000), (4000, 6500)} respectively. For these param- 
eters not more than 10% of the atoms is excited during 
each pulse, and thus the conditions of the adiabatic elim- 
ination arc fulfilled. The temperature was determined by 
fitting the calculated distribution of fermions to a ther- 
mal distribution. As one can observe, a final temperature 
T ~ 0.008TV may be reached within 8 s. The losses asso- 
ciated with the cooling process do not exceed 2% and are 
slightly larger for the laser-cooled component, due to the 
removal of the long- living excited atoms (inset of Fig.^l. 

We have analyzed the effect of the losses on the laser- 
cooled gas. To this aim we have considered an initial gas 
at T = 0.008Tf (end of the cooling process in Fig. [I}, an d 
compared two different cases: (i) the laser is turned off, 
and the gas is heated due to the background collisions, 
(ii) the laser is turned on, and the cooling pulses of the 
last stage are continuously applied. Fig. |2 presents the 
evolution of the temperature of both species for these 
two cases. As observed, the laser cooling compensates 




[s] 



FIG. 3: Heating of the degenerate two-component Fermi gas 
due to background collisions with a rate 7b g = 10 Hz. The 
plot shows the temperature of the two components obtained 
in the simulation (darker and lighter thin solid curves), the 
analytic result of Ref. fl9|l calculated for a homogeneous sys- 
tem (thick solid line), and the data calculated from Equation 
@ (thick dashed line) . 
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for the heating induced by the creation of holes in the 
degenerate distribution. Hence, it helps to maintain the 
degenerate gas for a relatively long time in the trap. We 
have verified that the continuous application of the laser 
cooling does not lead to substantially larger losses. 

Finally, we have analyzed the results for the heating 
due to background collisions with the help of an ana- 
lytical model. For the background losses that decrease 
the population of the trap levels in the way described 
by Eq. JHJ, one can calculate the evolution of the tem- 
perature analytically [l9| . To this end, it is sufficient to 
observe that the Eq. ||SJ do not change the mean en- 
ergy per particle e, which has to remain constant in the 
absence of cooling: e(N(t), T(t)) — const. At low tem- 
peratures T = T/Tp <C 1, the mean energy per particle 
in a harmonic trap is given by e(N,T) = (3/ A)Ep(N) + 
(TT 2 /2)(k B T) 2 /E F (N). Using £ (N(t),T{t)) = e(No,T ), 
where Nq is the initial number of atoms, and To is the 
initial temperature of the gas, we obtain 

f(f\- ! EF{Na) f 2 i 3 ( Ef{Nq) A en 

For a large system Ep(N) ~ hwUiN) 1 ^, and N(t) ~ 
iVo exp(— Tbgi), hence 

f (t) = ^/e^*/ 3 f 2 + 3 ( e W/3 - 1) /2tt 2 . (10) 

Fig. presents the evolution of the temperature for a 
rate 7b g = 10 Hz, which is much larger than the typ- 
ical experimental rates, and was employed to reduce 



the amount of time-consuming numerical calculations. 
Fig. compares our numerical results, the predictions 
of the analytical formula of derived for a homoge- 
neous system, and the data calculated from Eq. © with 
Ep(N) ~ hui(( l oN) 1 /' i — 2), which is more appropriate 
for a finite size system [2^. The analytic curve fits very 
well to the numerical data. The small discrepancy for 
larger temperatures originate from the fact, that the ex- 
pression for the mean energy, which we use, is valid up 
to order (T/Tp) 2 . From the plot, it is clearly seen that 
the increase of the temperature in a trapped gas is much 
slower that in the homogeneous case. 

In conclusion, we have studied the laser cooling of 
trapped two-component Fermi gases. Our method ex- 
ploits the dynamical adjusting of the effective sponta- 
neous emission rate in Raman cooling. Using Monte 
Carlo simulations, we have shown that laser cooling is 
able to cool the two-component Fermi system below 
0.01 Tp. We have also discussed the losses which may af- 
fect the laser-cooled gas. In this context, we have shown 
that our laser-cooling scheme can be employed to main- 
tain a two-component Fermi gas at a fixed temperature 
in the presence of background collisions, without signifi- 
cant additional atom losses. Finally we have derived an 
analytic formula for the temperature of a trapped Fermi 
gas heated by background collisions. We have shown that 
the heating rate is smaller than in the homogeneous case. 
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